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ABSTRACT 


Tupaz  et  al .  C^  9  78)  formulated  a  model  to  study  the  be- 
havior of  waves  in  an  unstable  jet  which  varied  downstream. 
The  linearized  barotropic  vorticity  equation  was  solved 
numerically.   Waves  were  forced  with  a  fixed  period  on  the 
eastern  boundary  and  a  radiation  condition  was  applied  on 
the  western  boundary. 

In  this  thes i s»  cycl i c  boundary  conditions  are  used  on 
the  eastern  and  western  boundaries.   The  numerical  solutions 
show  amplitude  growth  in  time  which  is  approximately  expo- 
nential.  The  solutions  are  normalized  and  the  wave  struc- 
ture is  obtained  during  three  time  segments.   The  wave 
structures  are  not  the  same  during  these  segments,  but  they 
are  very    similar  in  areas  of  large  amplitudes.   In  these 
regions  the  behavior  is  \/ery    similar  to  that  obtained  by 
Tupaz  et  al .  (1978) . 
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I .       INTRODUCTION 

A    jet,    regardless    of    its    latitudinal    position    in    the    at- 
mosphere,   may    contain    regions    of    large    vorticity    gradients 
where    the    necessary    condition    for    barotropic    instability    is 
locally    satisfied.      Synoptic-scale    moving    disturbances    occur- 
ring   at    the    level    of   a    moderately    strong    easterly    jet    south 
of   the    Tibetan    high    near    200    mb    level    appear    to    arise    from 
barotropic    instability    of    the    mean    flow    (Krishnamurti,    1971a, 
1971b).       If    the    observed    disturbances    are    the    result   of    baro- 
tropic   instability,    they    will    extract    energy    from    the    mean    flow. 

Many    investigators    have    studied    the    linear    stability    of 
barotropic    zonal    flows    over    the   years.      Kri s hnamur ti ' s    obser- 
vational   and    Colton's    (1973)    numerical     studies    both    showed 
that    the    zonal    variation    of    the    jet    apparently    had    significant 
effects    on    the    dynamic    behavior    of    the    transient    disturbances. 

Tupaz    et    al.    (1978)    further   extended    the    study    of    the 
zonal    variation    of   the    basic    flow.      A    numerical    model    was 
developed    based    on    the    linearized    non-divergent    barotropic 
vorticity    equation    on    a    beta-plane.      This    study    examined    the 
dynamics    of    transient    barotropic    waves    in    a    region    of    variable 
mean    wind.      The   mean    zonal    wind    was    an    easterly    hyperbolic 
secant-squared,    or    Bickley,    jet    and    the   mean    meridional    wind 
was    derived    so    that    the    mean    flow    would    be    non -di ve rgen t . 
Using    the    numerical    model,    a    periodic    forcing    was    applied    on 
the    inflow    (eastern)    boundary    to    generate    waves    which    were 
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allowed  to  propagate  through  the  channel  and  out  through  the 
western  boundary.   As  the  waves  moved,  they  reacted  to  the 
local  stability  properties  of  the  mean  flow,  whereas  at  each 
point  the  fields  varied  periodically.   The  results,  which 
were  obtained  from  long  term  integration,  were  compared  with 
the  parallel  flow  theory  by  constructing  a  simple  mechanistic 
analytical  model  which  incorporated  the  local  stability  con- 
cept of  the  parallel  flow  theory. 

In  this  thesis, the  model  developed  by  Tupaz  et  al.  (1978) 
is  modified  by  changing  the  boundary  conditions.   The  quanti- 
ties on  the  eastern  and  western  ends  of  the  channel  are  set 
equal  to  enforce  cyclic  continuity.   This  allows  a  wider 
variety  of  wind  profiles  to  be  considered  as  long  as  they 
are  periodic  in  x.   The  use  of  these  boundary  conditions  will 
also  improve  our  understanding  of  the  mathematical  behavior 
of  waves  in  a  mean  flow  with  downstream  variation. 
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II  .   BASIC  EQUATIONS 

A.   GENERAL  FORMULATION 

This  study  employs  the  barotropic  vorticity  equation  to 
describe  motion  in  a  periodic  channel.   The  development  which 
follows  Tupaz  et  al .  (1978)  begins  with  the  barotropic  vor- 
ticity equation  on  a  beta  (B)  plane: 


If  +  »  ¥t  +  •  If  +  8»  -  i 


DfC    , 


(2.1) 


where 


5    = 


_3v_ 
3x 


in 

3y 


(2.2) 


Here,     8      is    the    constant    value    of    the    north-south    gradient    of 
Coriolis    parameter    computed    at    10°    latitude.      Q      is    a    forcing 
function    which    represents    the    non- barotropi c    and    diabatic 
effects    which    are    required    to    maintain    the    basic    state    vor- 
ticity   field,    and      D£      is    a    frictional    coefficient.      The    mo- 
tion   is    assumed    to    be    non-divergent    so    that    a    stream    function 
(ijj)    can    be    introduced   which    is    defined    by 


dy    '  3x 


(2.3) 


The    vorticity    becomes 


5  ■     vij>    . 


2.4 
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The    non-divergent    barotropic    vorticity    equation    may    now    be 
written, 

&*-£*&♦& *£*  ♦»£■«-  v2*-     «.SJ 

The  linearization  of  eq.  (2.5)  is  accomplished  by  separat- 
ing the  stream  function  into  its  mean  (bar)  and  perturbation 
(prime)  quantities  of  the  form: 


i|>(x,y,t)  =  <Hx,y)  +  f  (x,y,t)  , 


(2.6) 


where       |  ip '  |    <    <    |  ij>  |     .      Substituting    (2.6)    into    (2.5)    and 
dropping    the    quadratic    terms    in      \p'       leads    to 


3      „2  .  ,  3^    8      n2  ,  ,     .     3ij)    3      _2  ,  ,  3'J;'     3      n2T 


3*'     3      v2-  3^'     . 

+  7x      "37   V    *         6  ~3~7 


DfV2<Jj 


(2.7) 


It  is  assumed  that  Q '  =  0  and  Q  is  defined  in  such  a  way 
that  ij)  is  a  steady -state  solution  of  (2.5).  Equation  (2.7 
can    also    be    rewritten    in    Jacob i an    form    such    that 


v2  ||'-  -  3n$.vV) 


where  a,b      =  -^—  ■*—•  -   -~—  tt— 

\-U  3x    3y         3y    3x 


3^ 


-    DfvV     , 


(2.8) 


Equation    (2.8)    is    a    Poisson    equation    for    the    tendency    of    the 

3iij  ' 
perturbation    stream    function    ( ■*■&    )  . 

6  T. 
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Equation  (2.7)  is  solved  in  a  channel  defined  by 


-  D  <  y  <  D 


xy  1  x  1  XE 


The  boundary  conditions  of  the  north-south  walls   y  =  +  D 
are 


*'  -  0  . 


(2.9) 


Cyclic    boundary    conditions    are    used    in    the    east-west    direc- 
tion.     Tupaz    et    al.    (1978)    applied    a    time-periodic    forcing 
on    the    eastern    boundary    to    introduce    waves    into    the    region 
from    the    east,    and    used    a    radiation    condition    for    both    ten- 
dency   and    vorticity    at    the   western    boundary    to    allow    the 
waves    to    move    smoothly    through    the    boundary. 

B.       ZONALLY    VARYING    BASIC    FLOW 

The    basic    velocity    field    is    defined    by    an    easterly    Bick- 
ley   jet 


u(x,y)    =    U(x)    sech2(^j)    -    UQ    = 


3y 


(2.10) 


Here,  d(x)   is  a  characteristic  length  scale  of  the  jet,  U(x) 
is  the  maximum  velocity  at   y  =  0,  and   U    is  a  constant 
velocity.   The  basic  flow  stream  function  is  specified  to  be 
constant  at  the  northern  and  southern  boundaries  (y  =»  +  D). 
Equation  (2.10)  can  be  integrated  to  give 


$(x,y)  =  U(x)  d(x){tanh(^y)  +  tanh( 


D 

dTTT 


)} 


UQ(y  +  D)  +  ip(-D)  , 


2.11) 
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where 


$(D)    +    2UQD   -    $(-D) 

u(x)  =  2dlD 


coth     (^y)     .         (2.12) 


Therefore,    if      d(x)      varies    slowly    in      x    ,    so    does    the    basic 
flow.      The    x-variation    for    the    characteristic    length    scale    is 

given    by 

(x_x    ) 
850    km    +    350    km{cos[2Tr      „    °    ]}     , 
d    =  L 

1200    km 


x      <    x    <    x 
o   -        -      o 

x      <    -    x„ 


x      >    x 
—     o 


(2.13 


Here,     x        is    the    longitude    where    the    x-variation    of    the    cosine 
o 

function    starts    and      L      is    the    wavelength    of    this    variation. 
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III.   FINITE  DIFFERENCE  EQUATIONS 

Equation  (2.7)  is  solved  by  finite  differences.  Centered 
time  differences  are  used  except  for  the  friction  term  which, 
is  approximated  with  a  forward  step.  The  Jacob ians  are  com- 
puted with  the  form  which  was  developed  by  Arakawa  (1966)  for 
non-linear  conservation. 

The  tendency  of  the  disturbance  stream  function  is  de- 
fined by 


Tt    _  r3^'t 
Ti,j  "  [Tt  ]i,j  ' 


(3.1) 


where    superscript      t      is    the    current    time    and    subscripts      i 
and      j      refer    to    the      x      and      y      grid    points,  respectively. 
In    the    following    finite    difference    equations,    At(l    hour)    is 
the    time    step.      Ax(375    km)    and    Ay(125    km)    are    the      x      and 
y      grid    point    intervals,    respectively.      Equation    (2.7)    is 
written    using    the    Leapfrog    finite    difference    scheme: 

v2ri,j =  -Ji-.j^2^  -j^^'.v2*] 


2Ax  V    *1  ,j 


,  , t+At  ,  , t-At  - 

iM     .  =    \b\      .  +    T.      .At 

l  'J  1  »J  1  .J 


(3.2) 
(3.3) 
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The  finite  difference  form  for  the  Laplacian  of  the  tendency 
is  given  by 


V2^ 


T*.,  .-2T*  .+T.  ,  .    T*  .x1-2T.  .+T.  .  , 


i  >j 


Ax 


Ay 


(3.4) 


and    for    the    Jacobian    terms    it    is    given    by: 


© 


jftj  =  4AW[(^i+i.j^i-i,j)UCj+r^,j-i)  <3-5) 


■  (*i.j+r*i.j-i)(«i£i.j"5i!i,j>]  • 


© 


Jl.J    =    4AxAy  [i^i+l  ,j(ci+l  ,j+l    "    ^i+l.j-^  (3,6 


^1-l,JUi-l,J+l    "    'l-l.j-l* 


*i  ,j+1  (ci+l  ,j+l    "    ci-1  ,j+l  5 


■  t 


,t 


+   *i,J-lUi+l,J-l    "    ^i-l,j-l)] 


and 


4AxAy  C?i  ,j+l  ^i+1  ,j+l 


,t 


'l-1.j+lJ 


*1-1.J-1> 


ci+i,j(*f+i,j+i  "  ^i+i ,o-i J 


+  ^-i5j(*i-i,j+i 


^i-i,j-i)]   ■ 


(3.7) 


1  7 


where 


i  (j 


CD 


i  ..i 


(3.3) 


i  ,,1 


3"  •  -[^'»v  $1  1S  expanded  in  the  same  manner  as  for  (3.8). 
This  Jacobian  conserves  both  kinetic  energy  and  enstrophy  in 
the  non-linear  formulations  (Arakawa,  1966). 

A  forward  time  step  is  used  for  the  initial  time  step  and 
every    30  time  steps  thereafter  to  damp  the  internal  gravity 
waves  produced  by  the  leapfrog  method.   A  Matsuno  finite  dif- 
ference scheme  is  tested  in  place  of  the  forward  time  step 
and  the  results  show  that  the  Matsuno  scheme  provides  further 
damping,  but  that  the  change  is  not  significant. 

The  boundary  conditions  employed  by  Tupaz  et  al .  (1978) 
relieved  the  problem  of  reflection  of  incident  waves  from  the 
boundary  back  into  the  interior  region  (Pearson,  1974).   An 
Eul er- backward  finite  difference  scheme  was  employed  every 
time  step  to  provide  additional  damping  of  the  short  waves. 

The  frictional  terms  are  evaluated  at  the  previous  time 
step  in  order  to  insure  linear  computation  stability  (Haltiner, 
1971). 

The  Poisson  Eq.  (3.1)  is  solved  for  the  tendency  using  the 
P0I5S  subroutine  which  is  a  cyclic  reduction  algorithm  for 
solving  block  tridiagonal  systems  of  arbitrary  dimensions  with 
Dirichlet  boundary  conditions  (Swarztrauber  and  Sweet,  1973). 

The  subroutine  POISS  solves  the  linear  system  of  equations 
of  the  f orm : 
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A.X.  ■,      .+B.X.  ,+C-X.^,  .+X.  .  ,-2X.  .+X.  ._,,  =  Y.  . 


for   I  =  1,2,  ...  M   and   J  =  1,2,  ...  N 


(3.9) 


The  east-west  boundary  conditions  for  all  quantities  are 
represented  by 


X-    =  X.  ,    i  =  0,  I 

1,0      1,1 


;i,J  +  l  =  Xi,l    J  =  °«  J  +  1 


(3.10) 


It    should    be    noted    that    Eq .    (3.9)    is    not    in    the    correct    form 
which    satisfies    the    model.      Therefore,    the   model    must    be 
flipped    prior    to    entering    the    solver    and    then    flipped    back 
upon    exiting   with    the    solved    eigenvalue-solutions. 
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IV.   LOCAL  STABILITY  OF  BASIC  FLOW 

A  parallel  flow  model  of  Williams  et  al.  (1971)  is  used 
to  determine  the  local  growth  rate  of  the  variable  mean  flow. 
This  is  taken  from  Tupaz  et  al.  (1978)  and  then  transposed 
to  fit  the  basic  flow  model  used  for  this  study  in  order  to 
gain  some  insight  on  the  stability  characteristics  of  the 
mean  flow.   This  model  is  hereafter  referenced  as  the  par- 
allel flow  model.   If   3ij;/3x   is  set  equal  to  zero  in  (2.5) 
the  governing  equation  of  this  model  is 


[Tt  +  °  Ix  +  Df]  V^'  +  [e 


a2u ,  a  if)1 

- 2"] 


3y 


3x 


0  , 


(4.1) 


u(y)  =  -U  sech2(^-)  +  UQ 


where   U   is  a  specified  constant  that  scales  the  magnitude 
of  the  central  velocity  of  the  Bickley  jet  (y=0).   The 
characteristic  length   d   is  also  a  specified  constant. 
Equation  (4.1)  is  solved  with  the  initial  value  technique 
used  by  Williams  et  al .  (1971).   Assuming  that  all  perturba- 
tions are  periodic  in  x,  Eq.  (4.1)  is  Fourier  transformed  in 
x  and  the  resulting  equations  are  integrated  until  the  solu- 
tion grows. with  a  constant  exponential  rate.   This  approach 
gives  the  growth  rate,  phase  speed  and  wave  number  structure 
of  the  most  unstable  meridional  mode  as  a  function  of  the 
zonal  wave  number  k  . 
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In  general,  these  equations  have  a  set  of  normal  mode 
solutions  and  a  continuous  spectrum  of  solutions  (Case,  1960; 
Pedlosky,  1964  and  Yanai  and  Nitta,  1968).   However,  only 
the  normal  mode  solutions  can  give  significant  growth  and 
the  most  unstable  mode  will  dominate  after  a  sufficient 
period  of  time. 

Figure  1  shows  the  growth  rate  (n)  corresponding  to  the 
most  unstable  wavelengths  (L)  as  a  function  of  x  based  on 
the  parallel  flow  model.   The  parameters  for  the  mean  flow 
are  given  in  Section  V.   The  largest  growth  rate  occurs  at 
x=0  where  the  jet  achieves  its  maximum  central  velocity 
[u  =  -  30  m  sec"  ].   In  the  inflow  and  outflow  region  for 
x  <  -  11,100  km  and  x  >  11,100  km,  respectively,  the  growth 
rate  becomes  negative  and  approaches  the  linear  frictional 
damping  rate.   The  most  unstable  wavelengths  range  from 
3650  km  at  x=0  to  4600  km  in  the  inflow  and  outflow  regions. 
This  variation  follows  the  variation  in  the  y-scale  of  jet 
d(x)  which  is  given  by  (2.16)  (Tupaz  et  al .  ,  1978). 

Ei gensol uti ons  of  the  most  unstable  discrete  mode  for 
two  jet  profiles  are  shown  in  Fig.  2.   Only  the  lower  half 
of  the  y  domain  is  shown  because  the  solutions  are  symmetric 
about  y=0.   Both  solutions  show  barotropic  growth  as  the 
phase  tilt  is  opposite  to  the  shear  in  u.   The  maximum  tilt 
occurs  at  x=0  where  the  growth  rate  is  the  largest.   In 
addition,  the  ei gensol uti on  amplitude  is  a  maximum  at  y=0  for 
x=0,  but  at  x  =  9314  km  the  maximum  has  shifted  to  y  =  +  800 
km . 
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V.   NUMERICAL  RESULTS  FOR  EXPERIMENT  I 

The  numerical  investigation  of  the  variable  mean  wind 
model  comprises  two  experiments.   Experiment  I  is  the  princi- 
pal one  and  it  employs  the  same  friction  that  was  used  by 

Tupaz  et  al.  (1978).   A  friction  coefficient  of   Df  =  0.15 

-  5    - 1 
x  10    sec    is  used  which  is  equivalent  to  an  e- folding 

decay  time  of  approximately  7.7  days.   The  second  experiment 

is  the  frictionless  case  and  will  be  discussed  in  Section  VI. 

The  following  parameters  are  used  in  Experiment  1. 


2D 

Uo 

At 


=    4000    km,   x       -   xw    =    40,125    km,    0(0,0)    =    -30    m    s 
=    0,    L    =    36,000    km,    Ax    =    375    km,    Ay    =    125    km    , 


1 


=    1    hr,   x 


+    18,000    km 


Figures    3    and    4    show    the    stream    function      $(x,y)      and    the 

zonal    velocity      u(x,y)      which    are    derived    from    the    above 

parameters.      The    model    has    108    grid    intervals    in    the    x-direc- 

tion    and    32    intervals    in    the    y-direction.      Tupaz    et    al.    (1978) 

used    the    same    domain,    but    his    mean    flow   was    not    the    same    at 

x    =XE      and     x    =xy      However    in    this    study,    the    scale    of 

variation    of   the    jet    is    reduced    to    insure    that    the    basic 

zonal    velocity   would    be    periodic    in    the    inflow    and    outflow 

regions.      Equation    (2.13)    specifies    that    the    jet    scale,    d(x), 

varies    between    500    km    at      u(x,0)  and    1200    km    at      u ( x ,  0 )     . 

max  'mm 

The  initial  disturbance  field  is  of  the  form: 
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,.,  *         sinTrycos2TTl0x 

♦  •(x.y.o)    -  2D       xr  -x„ 

E    W 


(5.1) 


A  wave  number  of  10  is  chosen  as  a  result  of  the  average  wave- 
length which  was  obtained  by  Tupaz  et  al .  (1978).   The  initial 
conditions  and  the  jet  flow  are  both  symmetric  about  y=0  and 
this  symmetry  is  maintained  throughout  the  integration. 
Therefore,  only  the  lower  half  of  the  y  plane  will  be  depicted 
in  the  results  of  both  experiments. 

The  integration  of  the  numerical  model  is  extended  for  a 
period  of  170  days.   A  close  examination  of  wave  structure, 
phase  tilt,  period,  wavelength,  and  spatial  growth  is  made 
for  the  segments  80-87  days,  125-133  days  and  150-157  days. 
Overall  growth  occurs  during  the  entire  period  and  is  expected 
since  the  waves  travel  through  the  region  and  then  reappear 
at  the  input  boundary.   It  was  expected  that  the  overall 
growth  would  be  exponential  in  behavior.   Figure  5  shows  the 
root  mean  square  for  the  entire  domain  of  \p'       plotted  on  a 
log  scale  vs  time.   If  the  growth  were  truly  exponential,  the 
curve  in  Fig.  5  would  be  linear.   But  in  fact  the  curve  is 
quasi -linear  with  departures  from  the  linear  behavior  occur- 
ring periodically  about  every    56  days.   Figure  5  shows  the 
first  departure  occurs  near  t  =  21  days  and  continues  for  28 
days.   It  should  also  be  noted  that  the  growth  rate  is  nega- 
tive from  t=38  days  to  t=47  days.   This  departure  from  the 
linear  appears  again  at  t  =  7 5  days  for  25  days  and  again  at 
t = 1 35  days  for  26  days. 
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Fig.   5.     Experiment   I 


'Jj  I  as   a   function  of  time,  where 

a  ve 


\\ii\         is   the  root  mean  square  of  iV 
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The   (fi1   field  at  t  =  30  days  is  shown  in  Fig.  6.   An  entire 
train  of  barotropic  waves  actually  exists  throughout  the  length 
of  the  channel,  but  the  waves  upstream  of  x  =  -3500  km  and  in 
the  outflow  region  are  not  s hown . because  of  their  relatively 
small  amplitude.   This  is  the  same  type  of  wave  behavior  ob- 
served  by  Tupaz  et  al.  (1978). 

The  phase  is  calculated  for  the  segments  mentioned  earlier. 
This  is  accomplished  by  observing  each  grid  point  in  the   i|»' 
field  and  recording  the  first  three  times  it  passes  through 
zero  in  either  direction.   A  seven  day  period  (P)  is  used  in 
order  to  insure  that  a  full  period  is  recorded.   Tupaz  et  al . 
(1978)  observed  a  periodicity  of  3.25  days.   The  actual  phase 
for  each  grid  point   is  computed  using  the  following  formula: 


360U(N-N  ) 
o  ' 


(5.2) 


where   H  represents  the  first  time  it  passed  through  zero  and 
N    is  the  starting  point  in  time  of  the  sign  calculations 
for  each  time  segment.   Although  the  wave  continues  to  grow 
during  the  time  segment,  the  period  should  not  change  appre- 
ciably.  Figure  7  represents  an  analyzed  section  of  the  phase 
tilt  for  the  segment  80-87  days.   We  note  that  all  the  waves 
shown  tilt  oppositely  to  the  mean  wind  shear,  which  is  neces- 
sary for  barotropic  instability.   The  maximum  tilt  is  observed 

downstream  of  x=0  where  the  jet  velocity  is  maximum. 

The  periods  at  y  =  0  and  y  =+750  km  for  all  three  time  seg- 
ments are  shown  in  Figs.  8-10.   Although  the  periodicity  is 
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not  constant,  it  is  about  77  hours  (3.2  days).   This  is  also 
representative  of  the  entire  channel.   Deviations  near  the 
inflow/outflow  boundaries  could  be  attributed  to  the  small 
amplitudes  in  this  region.   Figures  10a  and  b  show  an  increase 
in  period  at  x  =  -6500  km  to  93  hours  (3.87  days)  and  90  hours 
(3.75  days),  respectively,  and  a  decrease  at  x  =  6500  km  to 
39  hours  (1.65  days)  for  a  and  b,  respectively.   Figure  8 
shows  only  an  increase  at  x  =  6500  km  to  100  hours  (4.17  days) 
and  93  hours  (3.87  days)  for  a  and  b,  respectively.   This 
abnormal  behavior  could  be  connected  to  the  departure  of  the 
growth  rate  from  the  linear  behavior  illustrated  in  Fig.  5. 

Figure  11  is  an  overview  of  the  phase  tilt  for  all  three 
time  segments.   Although  they  are  not  the  same,  they  are 
similar  in  the  region  from  x  =  3875  to  x  =  -12000  km.   They 
tend  to  be  different  in  the  region  where  the  amplitude  is 
smal 1  . 

The  envelope  of  the  wave  packet  is  evaluated  at  y  =  -750 
km  where  the  wave  disturbance  amplitude  is  large.   The  nor- 
malized \p'       field  is  used  because  the  waves  are  still  grow- 
ing and  results  are  something  close  to  a  sinusoidal  variation. 
This  envelope,  <^'(x)>,   is  obtained  by  recording  the  average 

of  the  magnitudes  of  the  maximum  and  minimum   ti/ '  /  I  ii  I 

a  y    '  > r  *  a ve 

values  that  occur  at  each  longitude  over  a  period  of  7  days. 
Here,   j  ip  [      is  the  root  mean  square  \p'     .   Figure  12  shows 
the  wave  packet  of  all  three  time  segments  for  y  =  +  750  km. 
Note  that  the  maximum  amplitude  in  Fig.  12  is  near  the  point 
where  the  parallel  flow  growth  rate  drops  to  zero  in  Fig.  1. 


35 


V. 

% 

""         'v. 

• 

•n.. 

* 

^«. 

/               Z 

'"             V  * 

*                                   l        * 

1                                 -J          * 

\ •'"" 

*        * 

V 

• 

» 

»% 

- 

\l  • 

^» 

• 

• 
i  * 

r0 

•          / 

Y 

A 

•  •      J 

V 

j* 

■ 

(\ 

1  •* 

n 

m 

/•  * 

/•  N 

/  *       * 

/     •  -^  * 

^y*-*** 

f*  J 

*       • 

« 

•       £ 

• 

- 

1/ 

• 

m  £ 

* 
• 

(   / 

■ 
• 

«        / 

• 

[m 

• 

1 

J      . 

1  '   \        1 

i        i 

1 

en 

i 


I 


CT1 
I 


CxJ 

I 


I 


CO 

r— 
I 


I 


at  r-~. 

JZ   LT> 

■+->   r— 

1 

s-  o 

O   Uf> 

<+-   r— 

£.*— 

^L    O 

- - 

O 

IT)   "O 

r»-   c 

,  "S 

+  1^ 

ii     i 

• 

>>  i 

"O 

C 

T3     00 

>> 

o   ju 

T3 

II 

CO 

>,n 

^~ 

C       1 

<d  Ln 

a;  c-j 

3    r- 

+-> 

t 

0)  ^ 

-Q   -Cl 

*— — * 

*. ^ 

E 

-(-> 

.*: 

r—         *» 

»P"  - — ■» 

o 

+-> 

o 

• 

o 

0)      . 

1 — 

00 

— * 

ia 

-C    oo 

X 

a.  >> 

fO 

<<-  -a 

o 

i — 

2  CO 

OJ     1 

•r-    O 

>  CO 

s_ 

ai  ^ 

>  <a 

O'-' 

•    00 

i—i    4_) 

c 

-m  aj 

C    E 

<D   cn 

E    QJ 

•I—     00 

s_           • 

01  ai  co 

a-E   >, 

X  •(-    ro 

LU   +J  T3 

o 
in 


o 
o 


o 

LO 


o 

LT) 
I 


o 

o 


o 

LO 


(3Se^d)  0e 


36 


„     CM  i— 


en 


!»«. 

Ln 

QJ 

. — 

£ 

i 

+J 

o 

ID 

1. 

r— 

O 

4- 

* — * 

o 

i — 



a; 

T3 

-a 

o 

c: 

ts 

"3 

, 

^^ 

ra 

| 

u 

• 

•r— 

1 

s- 

• 

ai 



rrr 

3 

m 

C 

>> 

ra 

11 

-a 

+-> 

ro 

c"> 

4- 

, — 

O 

i 

»■ 

L[~> 

A 

CM 

.   o 


B 

V 

^ 

fl 

** 

c-> 

O) 

^^ 

I 

C3 

a. 

• 

o 

o 

• 

o 

. — 

• 

f— 

ai 

1 



> 

cz 

m 

1X3 

X 

ai 

<T3 

+-> 

-a 

0) 

.*: 

[•> 

u 

CO 

CTi 

<a 

i 

1 

Q.O 

CO 

o 

> 

• — -- 

<t> 

n3 

CM 

m 

■ ' 

r~ 

. 

LO 

y~* 

in 

4-> 

0)       • 

r-~ 

rz 

EO 

1 

O) 

cn 

S 

ai    1 

•r— 

in  — 

s_ 

CO 

ai 

ai    co 

l 

Q.  E    >> 

| 

X 

■r-     fO 

LU 

+->  T3 

CM 

I 


<(X)>> 


37 


Although  they  are  not  the  same,  it  should  be  noted  that  they 
are  downstream  of  the  max  jet  core. 

Tupaz  et  a  1 .  (1978)  derived  the  disturbance  kinetic 
energy  equation  to  better  understand  the  dynamics  of  baro- 
tropic  instability.   They  found  a  number  of  energy  terms  but 
observed  that  the  Reynolds  stress  term,  <(u'v'  3u/3y)>,  was 
the  only  dominate  source  for  disturbance  kinetic  energy.  This 
supported  the  large  amplitudes  and  spatial  gradients  of  i>' 
encountered  a  considerable  distance  downstream  of  the  jet. 
It  is  expected  that  this  energy  balance  will  also  hold  in  this 
study.   Since  this  is  a  region  of  tilt  opposite  to  the  shear, 
there  is  a  conversion  of  mean  flow  kinetic  energy  to  distur- 
bance flow  kinetic  energy.   Figures  1  and  2  illustrate  the 
large  energy  conversion  of  this  type  for  the  parallel  flow 
model  . 

Figures  13-16  show  the  phase  angle  0  of  the  wave  distur- 
bance as  a  function  of  y  for  three  different  time  segments. 
The  phase,  9  ,  for  the  most  unstable  wavelength  of  the  par- 
allel flow  model  is  also  included  (see  Section  4).   We  note 
that  the  waves  upstream  of  -10875  km  tilt  oppositely  to  the 
mean  wind  shear,  which  implies  barotropic  instability.   In 
one  case,  the  maximum  tilt  is  downstream  of  x=0,  which  is 
where  Tupaz  et  al .  (1978)  observed  the  maximum  tilt.   Between 
x  =  -10875  and  x  =  3750, the  phase  tilt  compares  closely  to 
the  parallel  flow  model  9  .  At  the  eastern  boundary  (Fig.  15b), 
9 2*  compares  rather  well  to  9  ,  but  9  *  and  9  *  show  exces- 
sive tilt  for  the  comparatively  weak  instability  in  this  re- 
gion.  This  excessive  fluctuation  may  partially  be  explained 
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by  the  small  amplitudes  which  are  present  in  this  region. 
Near  the  outflow  boundary  (Fig.  16)  the  tilt  is  reversed  in 
6_   and   8,.   This  indicates  dynamic  stability  or  a  flow  of 
energy  from  the  disturbance  back  to  the  mean  flow.   This  can 
also  be  detected  in  Figs,  lib  and  c.   In  almost  all  cases,  we 
observe  that  the  tilt  of  the  wave  disturbance  qualitatively 
adjusts  to  the  local  stability  of  the  mean  flow.   Since  the 
parallel  flow  model  can  only  solve  the  most  unstable  discrete 
mode,  dynamic  damping  is  not  indicated  in  Fig.  15  for  the 
outflow  region.   Note  that  the  phase  tilt  is  unstable  on  the 
average.   This  is  because  the  waves  show  a  net  growth  on 
moving  through  the  region. 

The  wavelengths  for  the  numerical  model  are  computed 
using  the  formul a  : 


L  = 


2tt  2Ax 

A8 


(5.3) 


where   A8   represents  the  change  in  phase  over   2Ax  . 
Equation  (5.3)  can  be  simplified  and  expressed  in  degrees 
resulting  in 


L  = 


360°  2Ax 

A6° 


(5.4) 


In  addition,  a  smoothing  process   is  applied  in  those  areas 
where  the  wavelength  varies  noticeably  to  give  a  more  accurate 
picture  of  the  wavelength  throughout  the  channel.   This  cal- 
culation is  made  by  choosing  a  questionable  point  (  O  )  in 
the  y  domain  of  the  phase  field  and  adding  180°  to  it.   This 
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new  phase  is  then  located  by  moving  in  both  x  directions.  The 
distance  between  the  new  phases  is  calculated  in  terms  of  the 
number  of  x  grid  points  and  multiplied  by   Ax   to  give  the 
corrected, or  smoothed, wavel engths .   Figures  17-19  show  the 
variation  of  wavelength,  L(x)  ,  for  all  three  time  segments 
at  y  =  _+  750  km.   L(x)    represents  the  smoothed  wavelength. 
Comparing  the  wavelengths  of  y=0  (not  shown)  and  y    =   ±   750  km, 
it  is  observed  that  the  wavelengths  are  quite  similar  for 
each  time  period.   The  maximum  wavelengths  are  within  +  4500 
km  of  x=0  and  tend  to  decrease  toward  the  boundaries  when 
the  wavelength  pattern  becomes  noisy.   This  is  the  same  region 
that  the  phase  tilt  (Figs.  15b,  16)  experiences  large  devia- 
tions from  the  parallel  flow  model  due  to  the  small  ampli- 
tudes.  Tupaz  et  al.  (1978)  observed  that  the  wavelength  at 
y  =  _+  750  km  was  larger  upstream  and  smaller  downstream  of 
x=0  than  at  y=0 . 

Figures  20-22  contain  the  spatial  growth,  m,  ., ,  computed 
from  the  numerical  model,  respectively.   Figure  23  is  the 
smoothed  version  of  m,  ^  and  the  spatial  growth  rate  (m)  of 
the  parallel  flow  model  obtained  by  using 


-n/C 


(5.4) 


where   C   is  the  phase  of  the  waves  and  n  is  the  local 


The  value  of  m,  .,  is  corn- 


growth  rate  (Tupaz  et  al.,  1978) 

puted  directly  from  the  wave  packet  envelope.   In  Fig.  22, 

* 
we  see  that  m ,  .,  has  a  larger  maximum  than  m  and  the  maximum 

is  shifted  slightly  upstream  from  the  jet  maximum.   All  three 
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curves  pass  through  zero  at  the  same  point  in  which  the  wave 
packet  envelopes  are  a  maximum  (see  Fig.  12).   Note  again  the 
noisy  pattern  near  the  eastern  boundary.   Tupaz  et  al  .  (1978) 
observed  these  same  features  except  that  the  m   curve  was 
skewed  to  the  left  with  respect  to  m  and  the  jet  maximum. 
This  was  the  result  of  the  tilt  structure  in  the  wave  lagging 
spatially  behind  the  expected  value  from  the  local  stability 
conditions.   The  (*)  represents  the  maximum  and  minimum  ob- 
served by  Tupaz  et  al.  (1978).   Note  the  similarity  of  the 

*      * 
maximum  and  minimum  to  m^   and  m,  „,  respectively.   The  reason 

* 
that  m.  .,  shows  more  damping  in  the  outflow  region  than  the 

parallel  flow  model  is  because  the  latter  can  only  give  the 

eigensolutions  of  the  most  unstable  mode. 
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VI.   NUMERICAL  RESULTS  FOR  EXPERIMENT  II 

Experiment  II  is  the  frictionless  case  and  employs 
the  same  parameters  and  calculations  procedures  that  are 
used  in  Experiment  I.   However,  although  the  integration  for 
this  experiment  is  extended  for  170  days,  only  the  time  seg- 
ments of  80-87  days  and  125-133  days  are  included  in  these 
resul ts . 

Figure  24  shows  the  root  mean  square   4> '   for  the  entire 
domain.  Note  that  the  total  growth  is  nine  orders  of  magni- 
tude greater  than  that  observed  in  Experiment  I  (see  Fig.  5). 
In  addition,  the  curve  in  Fig.  24  is  more  linear  in  appear- 
ance than  Fig.  5  but  still  experiences  the  departures  from 
the  linear  behavior  during  the  same  time  interval  that  were 
observed  when  friction  was  included. 

The  periods  at  y=0  and  y  =  +  750  km  show  very    little 
variation  and  closely  resembled  those  pictured  in  Figs.  7 
and  8.   The  predominating  period  which  is  representative  for 
the  entire  channel  is  the  same  as  Experiment  I  (77  hours  or 
3.2  days). 

Figure  25  is  an  overview  of  the  phase  tilt  for  the  two 
time  segments.   The  region  of  similarity  extends  from  x  = 
5000  km  to  x  =  -15,000  km.   When  the  amplitudes  are  small, 
the  behavior  becomes  erratic  as  observed  in  Experiment  I 
(see  Fig.  11). 
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The  wave  packet  envelope,  <^'(x)>,  is  shown  in  Fig.  2G. 
When  comparing  Fig.  2  6  to  Fig.  12,  it  is  observed  that  the 
curves  are  almost  identical.   Since  the  tilt  is  opposite  to 
the  shear,  there  is  a  conversion  of  mean  flow  kinetic  energy 
to  disturbance  flow  kinetic  energy. 

Figures  27-30  show  the  phase  angle  tilt  of  the  wave  dis- 
turbance y-structure  compared  to  the  most  unstable  wavelength 
of  the  parallel  flow  model  8  .   In  both  cases,  the  maximum 
tilt  was  at  x=0.   Comparing  them  to  Figs.  13-16  in  Experiment 
I,  we  see  a  similar  structure  except  for  Figs.  27  and  29a 
which  show  greater  tilt  between  500-750  km.   Note  again  the 
excessive  tilt  of  9  *  for  the  comparatively  weak  instability 
shown  in  Fig.  29b,  and  the  elimination  of  the  growth  of  0  * 
of  y  greater  than  -1000  km  in  Fig.  30  for  the  Figs.  15b  and 
16,  respectively,  in  Experiment  I.   The  curve  for  6  *  stresses 
the  dynamic  stability  which  can  also  be  observed  in  Fig.  25b. 

The  wavelength  patterns  for  y=0  and  y  =  +  750  are  wery 
similar  to  those  observed  in  Experiment  I.   The  maximum  wave- 
length occurred  within  +  4000  km  of  x=0. 

The  spatial  growth  rates  (Figs.  31-33)  are  near  mirror 
images  of  Figs.  20,  21  and  23.   The  growth  rates  are  a  little 
larger  than  in  Experiment  I  as  would  be  expected  because  there 
is  no  friction. 
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VII  .   SUMMARY  AND  CONCLUSIONS 

This  study  solved  the  model  developed  by  Tupaz  et  al  . 
(1978)  with  different  boundary  conditions.   They  examined 
the  behavior  of  linear  waves  in  a  zonally  varying  easterly 
jet  which  was  ba ro tropi cal ly  unstable  in  certain  longitudes. 
The  waves  were  forced  with  a  period   oj   on  the  eastern 
boundary  and  a  radiation  condition  was  applied  on  the  west- 
ern boundary  to  allow  the  waves  to  move  smoothly  out  of  the 
region.   After  a  period  of  adjustment  the  solution  at  every 
point  varied  with  period   oj  ,  but  as  each  wave  moved  through 
the  domain  it  experienced  growth  or  decay  in  response  to 
the  local  mean  fl'ow. 

In  the  present  study,  the  Tupaz  model  was  modified  by 
changing  the  boundary  conditions  to  cyclic  continuity  on  the 
eastern  and  western  ends  of  the  channel.   In  addition,  the 
mean  was  modified  slightly  so  that  it  satisfied  cyclic  con- 
tinuity.  Experiment  I  corresponded  to  Experiment  I  in 
Tupaz  et  al.  (1978)  and  Experiment  II  was  frictionless,  but 
otherwise  the  same  as  Experiment  I. 

In  Experiment  I,  the  root  mean  square  stream  function 
grew  with  time  because  as  each  wave  moved  through  the  channel, 
its  net  growth  was  added  when  it  reentered  the  region.   This 
growth  was  approximately  exponential  with  the  departures  from 
exponential  growth  having  a  period  of  about  56  days.   In  order 
to  compare  the  wave  behavior  in  this  study  with  that  of 


65 


Tupaz  et  al  .  (1978),  the  fields  were  normalized  to  remove 
the  overgrowth  in  time.   A  procedure  was  developed  to  deter- 
mine the  phase  and  amplitude  of  the  waves  during  a  particular 
time  segment.   The  following  time  segments  were  chosen  for 
analysis:   80-87  days,  125-133  days,  150-157  days.   The  wave 
structure  features  were  not  the  same  for  the  three  periods, 
but  the  features  were  very  similar  in  the  large  amplitude 
regions.   The  predominant  period  was  about  3.2  days  which  was 
very  close  to  the  forced  period  of  3.25  days  used  by  Tupaz 
et  al.  (1978).   The  largest  phase  tilt  was  found  near  the 
jet  maximum  and  in  two  of  the  segments  the  phase  tilt  reversed 
in  the  stable  regions  as  was  found  by  Tupaz  et  al.  (1978). 
The  amplitude  of  the  wave  envelope  was  a  maximum  near  the 
point  where  the  parallel  flow  growth  rate  was  zero.   The 
average  of  the  maximum  spatial  growth  rates  was  very    close 
to  the  value  obtained  by  Tupaz  et  al.  (1978).   This  verified 
their  result  that  the  spatial  variation  in  the  growth  rate 
could  augment  the  maximum  growth  rate.   The  downstream  wave- 
length was  somewhat  erratic,  but  it  was  generally  a  maximum 
near  the  jet  maximum  and  it  was  a  minimum  where  the  jet  was 
weak.   The  results  for  the  frictionless  case  were  similar, 
but  the  total  growth  was  much  larger. 

In  general,  the  wave  behavior  with  cyclic  boundary  condi- 
tions was  very    similar  to  that  obtained  by  Tupaz  et  al. 
(1978).   However,  the  solutions  in  this  thesis  were  more  dif- 
ficult to  analyze  because  the  wave  structure  in  the  three 
time  segments  was  not  the  same.   This  was  undoubtedly  related 
to  the  non-exponential  overall  growths.   Apparently  as  the 
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waves  came  out  of  the  barotropic  damping  region  and  entered 
the  jet  region  they  did  not  have  the  structure  to  begin  growth 
smoothly.   In  this  region  the  wave  phase  structure  was  par- 
ticularly erratic.   When  the  present  model  was  used  to  deter- 
mine wave  structure  and  behavior  with  other  wind  profiles, 
the  structure  should  have  been  computed  frequently  and  then 
the  typical  structure  could  have  been  obtained  by  averaging. 
This  study  should  be  extended  to  the  nonlinear  case,  by 
including  a  forcing  term  to  maintain  the  mean  flow  when  there 
are  no  waves  present.   When  waves  are  introduced  the  modifi- 
cation of  the  mean  flow  could  be  studied. 
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